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Sudden change of geometric quantum discord in finite temperature reservoirs 
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We investigate sudden change (SC) behaviors of the distance-based measures of geometric quantum discords 
(GQDs) for two non-interacting qubits subject to the two-sided and the one-sided thermal reservoirs. We found 
that the GQDs defined by different distances exhibit different SCs, and thus the SCs are the combined result of 
the chosen discord measure and the property of a state. We also found that the thermal reservoir may generate 
states having different orderings related to different GQDs. These inherent differences of the GQDs reveal that 
they are incompatible in characterizing quantum correlations both quantitatively and qualitatively. 
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I. INTRODUCTION 

Due to the emergence of quantum information science, the 
characterization and quantification of quantum correlations in 
a system have became one of the people’s research focuses 
in the past few decades Jll]. Nowadays, when we mention to 
quantum correlations, two prominent lines of research on this 
problem may immediately appear in our minds. The first one 
centers around the entanglement-separability paradigm, under 
which various forms of the entan^ment measures have been 
proposed and extensively studied Jj]. The second one is based 
on the noncommutativity of operators in quantum mechanics, 
and along this line people also presented a plenty of discord¬ 
like quantum correlation measures 0. 

It is well accepted that entanglement is responsible for the 
advantage of many quantum communication and computation 
tasks l3-0- It is also realized recently that quantum discord 
lH , which reveals quantum correlations from a different per¬ 
spective, plays a vital role in quantum protocols such as the 
deterministic quantum computation with one qubit 0, remote 
state prep aration 0 , quantum locking of classical correla¬ 
tions ilOl fTll] . quantum state broadcasting d, and quantum 
state merging lU^ . Moreover, the discord consumption 
has been linked to the quantum advantage for extracting in¬ 
formation via coherent interactions ifTsll , and meanwhile, it 
can also be connected to entanglement via some measurement 
processes d- 

Although both can serve as physical resources for quan¬ 
tum information processing, quantum discord and entangle¬ 
ment are in fact fundamentally different in many aspects. The 
most prominent one is that quantum discord may be increased 
by local operations 110 . while entanglement can only be in¬ 
creased by coherent operations. Moreover, when consider¬ 
ing their evolution for open quantum system, quantum discord 
and entanglement also exhibit distinct singular behaviors. For 
instance, quantum discord is more robust against decoherence 
than entanglement ifisl - l^ . it is immune to certain quantum 
noises during certain time intervals ll2ll - l^ . and therefore ex¬ 
hibits the frozen behavior which is impossible for entangle- 
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ment. Quantum discord exhibits sudden changes (SCs) due 
to the optimization procedure involved in its definition jldt - 
1^ . while entanglement undergoes sudden death which 
does not exist for quantum discord. It is also worthwhile to 
note that the sudden death of entanglement is independent of 
the chosen entanglement measures, as it occurs whenever the 
evolved state becomes separable. 

We investigate in this work the singular behaviors of three 
kinds of the geometric quantum discords (GQDs) for two non¬ 
interacting qubits subjecting to the independent thermal reser¬ 
voirs. We will show that for both the two-sided and the one¬ 
sided reservoirs, there are SCs being observed during the evo¬ 
lution of the GQDs, and the critical times for their occurrence 
are strongly dependent on the choice of the discord measures, 
which means that this phenomenon is the combined result of 
the chosen discord measure and the quantum state other than a 
property of the state itself. This statement is further confirmed 
by the relativity of different GQDs, i.e., different GQDs may 
impose different orderings of quantum states. 

We arrange this paper as follows. In SectionsHIlandlllll we 
recall the formula for the three GQDs, and present the analyt¬ 
ical solution for the master equation describing the evolution 
of the system. Then in Section ITVl we give a discussion of the 
SC phenomenon. Finally, Section[V]is devoted to a summary. 


II. DISTANCED-BASED MEASURES OF THE GQDS 

As we mentioned in SectionU the quantum discord in a sys¬ 
tem can be quantified from different perspectives. We adopt in 
this work three forms of the distance-based measures of GQDs 
||33-S21- They are defined via the minimal distance from p to 
Qqj i-e-j 

Da{p) = min da{p,x) (1) 

X6S7o 

where da designates the different distance measures one used 
in defining the GQDs. They are all well defined and can avoid 
the problem for the GQD based on the Frobenius norm S. 
We concentrate in the following on a bipartite system AB de¬ 
scribed by the density operator p, and the three GQDs are de¬ 
fined respectively via the trace distance 11341- the Hellinger 
distance |(3^|3^, and the Bures distance -|42|] . To facilitate 
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later description, we call them the TDD, HDD, and BDD for 
brevity. 

We first recall the corresponding measures for the trace dis¬ 
tance, the Hellinger distance, and the Bures distance between 
p and X, which are given respectively by 

di:{p,x) =11 P- X 111, 

dn{p,x) =\\ \fp-\/xh, (2) 

d^ip, x) = [ 2(1 - \/F{p,x))Y''^, 

where || • ||i and || • ||2 denote respectively the trace norm 
and the Frobenius norm, while F{p,x) = 
represents the Uhlmann fidelity. 

Then, the TDD, HDD, and BDD can be defined respectively 
via the minimal dT(p, x), '1 h(p, x), d^{p, x) 


Dt{p) = mm II p - X 111 , 

X 6 f 2 o 

Db_{p) =2 min \\ y/^- 

(3) 

T>b(p) = yJ{2 + V2) min(l- v/F(p,x)), 

where the set Hq is usually taken to be the classical-quantum 
state Ho = Efc(nfc ® lB)p{Tik ® Ib), and Hf, = Efe(nfc ® 
/B)v^(n^ (g) Ib), with {n^} being the projection-valued 
measurements. Moreover, the constants 2 and 2 -f i/2 be¬ 
fore min are introduced for the normalization of Dn (p) and 
Db{p) for the two-qubit maximally discordant states. 

The calculation of the TDD, HDD, and BDD is a hard task 
for general p, and analytical results exist only for certain spe¬ 
cial classes of states. First, for the two-qubit X state p^ which 
only contains nonzero elements along the main diagonal and 
anti-diagonal, the TDD is given by llj^ 


L>t(p ) = 


I 

Cmax - Cmin + C? ” ’ 


(4) 


where ^ 1,2 = 2(|p23| ± IPuI), 6 = 1- 2(p22 + P 33 ), Cmax = 
-f x\^}, and = min{C^,^|}, with xa 3 = 
2(pii + P 22 ) — 1- 

Second, for the 2 x n dimensional state p, and the decom¬ 
posed y/p = CtjXi 0 Yj, with {Xi : * = 0,1, 2, 3} and 
{Yj : j = 0,1,..., — 1} constituting the orthonormal op¬ 

erator bases for the Hilbert spaces TLa and TLb, the HDD can 
be calculated as ll^ 


where ||r ||2 = ^tid Zmax represents the largest eigen¬ 

value of the matrix with Z = (c^ )i^i_ 2 , 3 ;j=o,i, .. 

Finally, although there is no analytic solution, the calcula¬ 
tion for the maximum of F{p, x) can be simplified as lii^l 


6nax(p, x) = o max [ 1 - trA -f 2 V' Afe(A) 


(6) 


fc=i 


where Afc(A) denote the eigenvalues of A = yfp{u-5®lB)yfp 
in non-increasing order, with u being a unit vector in K^, and 
(T = {a^, a'^, a^) is the vector of the Pauli operators. 

From the above equations, one can note that even for the 
simple X state, there are optimization procedures involved for 
obtaining the GQDs, and this may induce SC behaviors of 
Dct{p). We will discuss them explicitly in the following text. 


III. SOLUTIONS OF THE MODEL 


The central system we considered consists of two qubits 
(labeled as S' = A, B) with large enough spatial distance, 
and the direct interaction between them can be ignored. We 
will discuss two different cases: (i) both qubits A and B are 
embedded in their own independent thermal reservoirs, and 
(ii) only qubit A (or B) is embedded in the thermal reservoir, 
while the other one is free of noise. 

For the qubit S subject to the thermal reservoir, the evolu¬ 
tion of p^{t) is governed by the master equation lii^l 

fc=l 


under the Markovian approximation. Here, 75 is the strength 
of the damping rate, {•, •} denotes the anticommutator, while 
Cf = y'n + lag and Cf = af being the rais¬ 

ing and the lowering operators) describe, respectively, the de¬ 
cay and excitation processes of the qubit S, with rates depend¬ 
ing on the temperature which is proportional to the average 
thermal photons n in the reservoir. 

The master equation (|2ll can be solved analytically. For con¬ 
venience of later presentation, we define 


9i = 


; + (n + l)p| 
2n + l 




b{1-Ps) 
2,71 -|-1 


( 8 ) 


Db{p) = 2(1 - ||r ||2 - Zmax), 


(5) then one can obtain that p^{t) takes the following form 


J 


dl PlliO) + Q 2 Pooi^) PSPloiO) 

^ V PsPoi(O) l-9fPn(0)-<7|Poo(0) 


(9) 


where = (*|p'^|j) in the standard basis {| 1 ), | 0 )}expanded 


by the eigenvectors of the Pauli operator ag, and the time- 
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dependent factor PS = Moreover, when n = 

0, the reservoir is at zero temperature, and the solution of Eq. 
(|9]l reduces to that given in Ref. . 

From the analytical expression for the single-qubit reduced 
density matrix (t), one can obtain the two-qubit density ma¬ 
trix p{t) for arbitrary initial state p(0) by using the procedure 
presented in Ref. 11^ . The diagonal elements are given by 

Pii{t) = gfgf pii(O) -f gfgf P22(0) -f P33(0) 

+ ^2 P44(0), 

P 22 {t) = gi^[(l - gf )pii(0) + (1 - gf )P22(0)] 

+ 9^[(i - gf )P33(o) + (1 - gf )p44(o)], (10) 

P33{t) = (1 - gf )[gf pii(O) -f gf P22(0)] 

+ (1 - <l 2 )kl P33i0) + gf P44(0)], 

P44(i) = 1 - - P22{t) - P33{t), 

while the nondiagonal elements are given by 

Pl2{t) = gfpBPl2(0) + gfpBP34(0), 

Pl3{t) = gfpAPl3(0) + gfpAP24(0), 

Pl4{t) = PAPBPuiO), P23{t) = PAPBP23i.O), (H) 

P24(t) =pa( 1 - gf)pi3(0) -f pa( 1 - gf)p24(0), 

P3i(t) = (1 - gf)pBPi2(0) + (1 - gf )pbP34(0), 

and the other nondiagonal elements can be written directly by 
using the Hermitian condition Pij(t) = p*i(t). 

For the special case of pa = Pb = P, we obtain the two- 
sided identical reservoir £ab = £a® £b, while for pA = p 
and Pb = 1 iPA = 1 and pB = p), it reduces to the one-sided 
reservoir £a (£b)- 

IV. SC OF THE GQDS 

Based on the solutions in Eqs. ( fTOb and (fTTT i for the system- 
reservoir coupling mode presented in Section [Till we begin 
to discuss decay dynamics of the three GQDs. We will show 
that they exhibit distinct singular behaviors, which include the 
completely different SCs and their relativity on characterizing 
quantum correlations. To be explicit, we consider initial two- 
qubit state of the following form 

IT-) = a|ll) + /3|00), (12) 

where a € [0,1], and = y/l — a^. The analytical expres¬ 
sions of p{t) for both the two-sided identical reservoirs £ab 
and the one-sided reservoir £a (or £b) can be written directly 
from Eqs. (fTOl i and (fTTI) . which are of the X form. 

We first consider the two-sided identical reservoir £ab- 
Fig. [T]is an exemplified plot of the '-yt dependence of Dt^{p), 
Db{p), and Db{p) for the initial state |T') with n = 0.6 and 
different values of . For this case, as P 23 {t) = 0, the TDD 
can be obtained analytically as 


(13) 



it 


FIG. 1: (Color online) yt dependence of Dt{p) (black), Db{p) 
(red), and Dh{p) (blue) for the initial state I'F) subject to the two- 
sided reservoir £ab with n = 0.6. The other parameters are given 
by = 0.3 (a), 0.5 (b), and 0.7 (c). The hollow circles denote the 
SC points. 


therefore it is symmetric with respect to = 0.5, and decays 
smoothly and monotonously with increasing yt for any . 

The BDD and the HDD are no longer the symmetric func¬ 
tions of = 0.5. As displayed in Fig. [T] while both 
Z?b(p) and I?h(p) still decay monotonously with increasing 
yt, there are also SCs being observed, i.e., they are nons¬ 
mooth functions of yt. In particular, the critical time ytc 
for the SCs and the times of SCs are determined strongly 
by the chosen discord measure and the form of the initial 
state. For the chosen parameters in Fig. [T] Db{p) decays 
smoothly for d^ = 0.3, while Dnip) experiences double SCs 
at ytc — 0.0115 and 0.287, respectively. For d^ = 0.5, the 
single SC of Db{p) {ytc — 0.478) occurs earlier than that 
of Db{p) {ytc — 0.5115), while for c? = 0.7, both Db{p) 
and 7 ?h(p) exhibit double SCs, where the first one of Db{p) 
{ytc — 0.0404) occurs shortly after Db{p) {ytc — 0.0066), 
and the second one of Db{p) {ytc — 0.5985) turns out to be a 
little bit earlier than that of Db{p) {ytc — 0.6113). 

Next we turn to discuss the cases of the one-sided reservoirs 
£a and £b. We displayed in Fig. |2]the yt dependence of the 
GQDs by the solid and the dash-dotted lines, respectively. As 
for these two cases, only one of the two qubits is exposed to 


D'y{p) = 2p^a \/1 — d^, 
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FIG. 2: (Color online) 'yt dependence of Dt(p) (black), Db(p) 
(red), and Dh(p) (blue) for the initial state |4') subject to the one¬ 
sided reservoir £a (solid) or £b (dash-dotted) with n = 0.6. The 
other parameters are o? = 0.3 (a), 0.5 (b), and 0.7 (c). The hol¬ 
low circles and squares denote the SC points. Moreover, the lines of 
Dt(p), as well as the lines of Db{p) with o? — 0.3, are overlapped 
for the £a and £b cases. 

the reservoir, the decay of the GQDs is slower than that for the 
two-sided reservoir case. First, for both 8a and Eb, the TDD 
is always given by 

Dt{p) = 2pa \/1 — (14) 

which is still a symmetric function about = 0.5, and de¬ 
cays smoothly and monotonously with the increasing 'yt. 

The BDD and the HDD may do not behave as smooth func¬ 
tions of jt. For the one-sided reservoir £a, they exhibit qual¬ 
itatively the same behaviors as those for the Eab case, and 
the only difference is that the critical times for the SCs are all 
obviously delayed. For c? = 0.3, the double SCs for Db{p) 
occur at 'ytc — 0.024 and 0.536, respectively. For = 0.5, 
the single SC for Db{p) and Db{p) occurs respectively at 
'ytc — 0.802 and 0.816. Finally, for = 0.7, the double SCs 
for Db{p) occur at 'ytc — 0.0611 and 0.966, while for Db{p) 
they occur at 'ytc — 0.01045 and 0.914. 

For the one-sided reservoir Eb, although the evolved den¬ 
sity matrix differs only in p 22 (f) P 33 [t) from the Ea case, 

the BDD and HDD are not exactly the same (cf. the solid and 
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FIG. 3: (Color online) 'yt dependence of Dt{p) (black), Db{p) 
(red), and Db{p) (blue) for the initial state I'T) subject to the two- 
sided reservoirs Eab (a), one-sided reservoir Ea (b), and £b (c), all 
with the parameters n = 0 and c? — 0.5. The hollow circles denote 
the SC points. 


the dash-dotted lines in Fig. |2|. For = 0.3, Db{p) ex¬ 
hibits completely the same 'yt dependence compared with that 
of the Ea case, while the double SCs of Db{p) {"ftc — 0.022 
and 0.652) are slightly different. Moreover, the single SC for 
both Db{p) and Db{p) disappears for = 0.5. Finally, for 
q;^ = 0.7, the first SC of Db{p) (yytc — 0.1305) and Db{p) 
{'ytc — 0.027) occurs later than that of the Ea case, while 
their second SC {'ytc — 0.822 and 0.724, respectively) occurs 
earlier than that of the Ea case. 

The decay rates of the GQDs for Ea and Eb may also be 
different. As showed in Fig. for = 0.5, Db{p) and 
Db (p) for the Eb case decay faster than those for the £a case 
in the whole ■yt region. For = 0.3 and 0.7, however, Ea 
and Eb give different decay rates of Db{p) and Db.{p) only 
during limited 'yt regions. For = 0.3, Db{p) for the Eb 
case decays slightly faster than that for the Ea case when 'yt S 
[0.022,0.652], while for = 0.7, Db{p) and Db{p) for 
the Eb case decay slower than those for the Ea case when 
'yt £ [0.0611,0.966] and 7 f £ [0.01045,0.914], respectively. 

We discussed in the above evolution of the three GQDs, and 
observed distinct singular behaviors such as the single and the 
double SCs caused exclusively by the reservoir. We discuss 
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in the following two limiting cases, i.e., the zero temperature 
(n = 0) and the infinite temperature (n —>■ oo) cases. 

For the zero temperature case, £1 = 0’ the only sponta¬ 
neous decay term leads to a purely dissipative process. In the 
long-time limit, this process drives the corresponding qubit 
to its ground state |1), thus the GQDs disappear for arbitrary 
initial state. In Fig. [3 we showed the ■jt dependence of the 
three GQDs for with = 0.5 and n = 0. For Sab, by 
comparing with Fig. [Tfb), one can see that they exhibit very 
similar behaviors, except that the decay rates are obviously 
decreased, and the critical times for the SCs ( 7 ^ — 0.485 and 
0.609, respectively) are also slightly delayed. For 6a, there 
is also single SC for both Db{p) (jic — 0.693) and Db{p) 
{•ytc — 0.609), which occurs earlier than those for the finite 
temperature reservoir (cf. Fig. [3b) and Fig. I3b)). Moreover, 
as showed by the blue line in Fig. Ob), the HDD is increased 
with jt during the region jt G [0.482,0.609] (which is a re¬ 
flection of the fact that the quantum discord can increase un¬ 
der the local operation on one party of the system ll4^ 1. while 
the TDD and the BDD are always decreased. This means that 
different discord measures may lead to different orderings of 
quantum states, and it confirmed again that what the discord 
reveals is in fact the combined result of the chosen discord and 
the quantum state other than a property of the state itself. Of 
course, the increase for HDD is slight and transient, and after 
'yt > 0.609, it decays to zero gradually. Finally, for 6b, the 
three GQDs still show qualitatively the same dependence 
with those for the finite temperature reservoirs, with however 
the decay rates are evidently decreased. 

For the infinite temperature case, we have £f = £2 , the 
decay and excitation processes occur at exactly the same rate, 
and the noise induced by the transitions between the two lev¬ 
els brings an arbitrary initial state into the maximally mixed 
one. For concise of the paper, we do not present the plots 
here. But the numerical results showed that for the initial state 
jT*) with = 0.5, all the three GQDs decay smoothly and 
monotonously with the increasing jot ( 70 = nj), and there 
are no SCs being observed for them. 

Finally, as the SCs displayed in Figs. (Tj [2] and |3 are ob¬ 
tained via numerical methods, one may wonder whether they 
are the real SCs or not, this is because sometimes it is possible 
that what one observes as a SC might be the result of a quick 
change that is actually not sudden when analyzed for smaller 
time intervals ifsotl . For the TDD, as its analytical expressions 
are given in Eqs. ( [T3] l and (O, it is evident that it does not ex¬ 
perience SC for the model considered here. For the HDD and 
BDD, as the square root of the density operator p cannot be 
derived analytically, analytical solutions of Db{p) and Db{p) 
cannot be obtained. But the changes observed in the three fig¬ 
ures are actually exactly sudden, and they are caused by the 
discontinuity of the optimal angle related to measurement op¬ 


erators {n^}. In fact, by writing 11^2 = ± u ■ it)/2, and 

the unit vector u = (sin 0 cos (j), sin 9 sin (j), cos 9), we found 
that for the considered state in this paper, both Db{p) and 
Db (p) are independent of the angle (j), but for the case of the 
single SC, the optimal 9 for both Db{p) and Db{p) is given 
by 7r/4 before the SC point, and it changes abruptly to 7r/2 
after the SC point. Moreover, for the case of double SCs, the 
optimal 9 changes abruptly from 7r/2 to 7r/4 and then back 
to 7r/2 at the SC points. All these correspond to SCs of the 
optimal {n^}, and therefore the changes observed for both 
Db (p) and Db (p) are actually sudden. 


V. SUMMARY 

In summary, we have investigated time evolution and the 
accompanying singular behaviors of the TDD, BDD, and 
HDD. To focus exclusively on the singular behaviors of them 
as caused solely by the thermal reservoir, the two qubits of 
the central system are assumed to be spatially separated far 
away from each other, and thus there are no direct interac¬ 
tions between them, that is, every qubit interacts with its own 
independent reservoir. 

By solving analytically the master equation describing the 
evolution of the two qubits, we analyzed dynamics of the three 
GQDs, and found that they are incompatible in characterizing 
quantum correlations, although they are all well defined from 
a geometric perspective. Our findings are illustrated through 
two distinct behaviors of the three GQDs. First, we found that 
the three GQDs may exhibit completely different SCs for both 
the two-sided and the one-sided reservoirs. The critical times 
for the SCs and the times of SCs are strongly dependent on the 
choice of the GQD measure and the form of the initial state. 
Different GQD measures have different SCs, and thus the SCs 
are in fact the combined result of the chosen GQD measure 
and the quantum state, but not the intrinsic property of a state 
itself. This is fundamentally different from the sudden death 
of entanglement, which is independent of the entanglement 
measure. Moreover, we also revealed the relativity of differ¬ 
ent GQDs. To be explicit, we found that the thermal reservoir 
may lead to a generation of quantum states manifesting differ¬ 
ent orderings, and this implies that different GQDs are incom¬ 
parable as their behaviors may not only be quantitatively but 
also be qualitatively different. 
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